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A comprehensive field theory is developed for superconductors with quenched disorder. We 
first show that the matrix field theory, used previously to describe a disordered Fermi liquid and 
a disordered itinerant ferromagnet, also has a saddle-point solution that describes a disordered 
superconductor. A general gap equation is obtained. We then expand about the saddle point to 
Gaussian order to explicitly obtain the physical correlation functions. The ultrasonic attenuation, 
number density susceptibility, spin density susceptibility and the electrical conductivity are used 
as examples. Results in the clean limit and in the disordered case are discussed respectively. This 
formalism is expected to be a powerful tool to study the quantum phase transitions between the 
normal metal state and the superconductor state. 
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I. INTRODUCTION 

The description of disordered many-electron systems is 
a difficult problem in modern theoretical physics. Lan- 
dau Fermi-liquid theory^ and then many-body perturba- 
tion theory^^'^ were introduced to deal with this problem. 
Considerable progress has been made within the frame- 
work of the latter. On the other hand, field-theoretic 
method has also been applied to the many-electron prob- 
lem, which has certain advantages over the traditional 
technique. It is relatively easy to include the quenched 
disorder in the formalism with the help of replica trick. 
More importantly, it is the natural language to describe 
any classical or quantum phase transition. It allows for a 
straightforward application of the renormalization group, 
implementing an old program of describing the various 
phases of many-body systems in terms of stable RG fixed 
points.* So far this program has been carried out for clean 
and disordered Fermi liquids, as well as disordered ferro- 
magnetic metals. ^'^ 

In the present paper we will develop a comprehen- 
sive field-theoretical method, or matrix field theory,® for 
gapped disordered spin-singlet superconductors. The for- 
malism is expected to be able to describe the quantum 
phase transition between the normal metal state and 
the superconductor state. Some similar techniques have 
been applied to describing spin-triplet, even-parity su- 
perconductors, but explicit quantitative expressions for 
the Gaussian propagators could not be given in the pre- 
vious paper and thus that description cannot be said to 
be complete.^ Here we will completely determine all the 
soft, or gapless correlation functions for the S = 0, spin- 
singlet case and obtain the corresponding transport prop- 
erties. The method can be generalized to evaluate other 
physical systems, like spin-triplet superconductors. Our 



results for the spin-singlet case coincide with earlier ones 
obtained by conventional methods. Our field theoretic 
methods, however, have the advantage that they can be 
easily generalized to describe quantum phase transitions. 
For example, in future publications we will use these 
results to (1) describe a metal-superconductor transi- 
tion in a dirty metal, without integrating out the soft 
fermionic degrees of freedom^ and (2) consider the same 
superconductor-metal transition from the superconduct- 
ing side of the transition.^ This second problem is non- 
trivial and interesting due to the numerous dangerous 
irrelevant variables at this phase transition. 

The paper is organized as follows. In Sec. II we give a 
field-theoretic formulation of the problem. In Sec. Ill we 
develop the theory to construct a saddle-point solution, 
obtain the gap equation of superconductivity and expand 
to the second or Gaussian order about the saddle point. 
In Sec. IV we show how to calculate the ultrasonic atten- 
uation coefficient, the number and spin density suscep- 
tibilities, and the electrical conductivity for both clean 
and disordered superconductors. In Sec. V we conclude 
with a general discussion of our results. In Appendix we 
give some technical points that are used in the paper. 



II. MATRIX FIELD THEORY 

A. Grassmannian field theory 

In general, a system of interacting, disordered electrons 
will be considered. The partition function of the system 
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Here the ip and ip arc Grassmann valued fields, and S is 
the action including three parts: 



S — Sq iSint + 'S'dis ) 



(2.2a) 



Where, with a (d + l)-vector notation x = (x, t) and 

J dx = Jy dx Jq dr, Sq describes free electrons with 
chemical potential /x, 

So = J dx^ tp„{x) (^~^-r + ^+l^ i>a{x) , 

(2.2b) 

S'int describes a spin-independent two-electron interac- 
tion, 

'S'int = - ^ j dxi dx2 ^ v{y.i - xa) 

Cl ,(72 

X Viai {Xi) V;<^2 {X2) V'<72 {X2) V'ai (sJl) (2.2c) 

and S'dis describes a static random potential ■u(x) cou- 
pling to the electronic number density, 

S'dis = - dx^^ u{x.)if;cr{x)ipcr{x) ■ (2. 2d) 

a- 

For further calculation, we assume the random poten- 
tial w(x) in Eq. (2. 2d) has variance. 



{«(x)«(y)}dis = 



1 



and is Gaussian distributed. 
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(2.3a) 



{•••}dis = j D[u]P[u]{...) , (2.3b) 

Np is the density of states at the Fermi level, and Te is 
the elastic scattering time. The disorder is quenched, so 
the replica trick^^ is used. With 



\nZ= lim (Z™- l)/m , 

m— >0 



(2.4) 



we consider, 
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jlD[r,r]^Ms] , (2.5) 
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where the corresponding action S equals to 
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It is also useful to get a Fourier representation with 

wave vectors k and fermionic Matsubara frequencies 
LOn = 27rT(n +1/2) by the following transformations: 
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r{x) = VT dr e^'^-^ i>„{x) , 
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dT e-'""'' ipa{x) , (2.7a) 



and 



dx e "^"^ tpna{-X.) , 



tpnaO^) = ^= I dx e"'-^Vn<T(x) 

vV 
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The procedure used here is similar to the one used in Ref. 
6, and we refer the reader to it for further details. 



B. Composite variables: Q-matrix 

Now we integrate out the Grassmann fields and rewrite 
the theory in terms of complex-number fields. As a 
first step, the resulting model can then be approximately 
solved by using saddle-point techniques. Later fluctua- 
tions about the saddle point will be considered. First we 
introduce a matrix of bilinear products of the fermion 
fields, 



I -V'itV:2t -^ntn -V'iTV'2i 

-V'u'^2t -tntn -V'u'02i 

V:^uV:^2T "01X^^21 '0uV'2i 

\ -V'itV'2T -V'iTV'2i -V'iTV'2i 
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■011^21 / 
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where all fields are understood to be taken at position 
X, and 1 = (ni,ai) with m denoting a Matsubara fre- 
quency and a a replica index, etc. The matrix elements of 
B commute with one another, and are therefore isomor- 
phic to classical or complex number-valued fields that we 
denote by Q. We use the notation a = 6 for "a is isomor- 
phic to 6" . This isomorphism maps the adjoint operation 
on products of fermion fields, which is denoted above by 
an overbar, onto the complex conjugation of the classical 
fields. We use the isomorphism to constrain B to the 
classical field Q by means of a functional 5 function, and 
exactly rewrite the partition function^ 

Z = j D[i,,i,] e^[^'*l Jd[Q]S[Q-B] 

' D[i;,^] e^'I'^'^^l / D[Q]D[A] e^I^^"^)! 



(2.9) 



= I D[Q]D[A] 



,A[Q,A] 



Here A is an auxiliary bosonic matrix field that plays 
the role of a Lagrange multiplier, and integrates out the 
fermion fields. 

It is useful to expand the 4x4 matrix in Eq. (2.8) in 
a spin- quaternion basis. 
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and analogously for A. Here tq = sq = I2 is the 2x2 unit 
matrix, and tj = —Sj ~ ^io'j, [j = lj2,3), with (Ti.2,3 
the Pauli matrices. In this basis, i = and i = 1,2,3 de- 
scribe the spin singlet and the spin triplet, respectively. 
An explicit calculation reveals that r = 0, 3 corresponds 
to the particle-hole channel (i.e., products tpip), while 
r = 1, 2 describes the particle-particle channel (i.e., prod- 
ucts ijj-ip or ijjijj). From the structure of Eq. (2.8) one ob- 
tains the following formal symmetry properties of the Q 
matrices,^ 



rQl2 ~ i ) 



(r = 0,3) , 



(-r+';Q2i > (^ = 0,3; i = 1,2,3) 



(2.11a) 

r^l2 — \^J r^21 ' V — "^i - — ■^t^^^J ) (2.11b) 
rQl2=rQ21 > = 1, 2) , (2.11c) 

;Qi2 = , (r = 1, 2; i = 1, 2, 3) , (2.11d) 



■ni — 1,-712 — 1 



(2.11e) 



Here the star in Eq. (2. lie) denotes complex conjugation. 

Now by using the delta constraint in Eq. (2.9) to 
rewrite all terms that are quartic in the fermion field 
in terms of Q, we can achieve an integrand that is bilin- 
ear in V' and ^jj. The Grassmannian integral can then be 
performed exactly, and we obtain for the effective action 
A 

A[Q, A] = AM + ^dis[g] + ^ Tr In (Gq ^ - it) 

+ dx tr (a(x) Q(x)) . (2.12) 

Here Tr denotes a trace over all degrees of freedom, in- 
cluding the continuous position variable, while tr is a 
trace over all those discrete indices that are not explic- 
itly shown. And 
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is the inverse free electron Green operator, with dr and 
c?x derivatives with respect to imaginary time and posi- 
tion, respectively, m is the electron mass, and ji is the 
chemical potential. We can see from the structure of 
the Tr In- term in Eq. (2.12) that the physical meaning 
of the auxiliary field A is that of a self-energy. The 
electron-electron interaction A^t is conveniently decom- 
posed into four pieces that describe the interaction in the 
particle-hole and particle-particle spin-singlet and spin- 
triplet channels.^ For the purposes of the present paper, 
we need only the particle-particle spin-singlet channel in- 
teraction explicitly to describe superconductivity. Sim- 
ilar to the BCS model we ignore the normal Coulomb 
repulsion in the particle-hole channels, and we also ig- 
nore the possibility of triplet superconductivity.^^ Then 



Ant[Q]=A„t 



-^int 



r— 1,2 ni ,712 ,m a 
X [tr {{rr^So)QZ-n,+mi^))] 

X [tr ((r.®so)g^:,,„,+„(x))] , (2.14) 

with r(°) the particle-particle spin-singlet channel inter- 
action amplitude, with F^'^^ < leading to superconduc- 
tivity. For the disorder part of the effective action one 
finds^^ 

^dis[Q] = j rfx tr (Q(x))' . (2.15) 

We will focus on the matrix elements qQ and iQ in 
disordered superconductivity states. From Eqs. (2.8) and 
(2.10) we find 

0Q12W = I [-V'iT(x)Vi2T(x) - Vu(x)Vi2i(x) 

+ Vii^(x)V'2i(x)+ViiT(x)V2T(x)] , (2.16a) 



iQui^) = -J- [-^iT(x)V2i(x) + Vu(x)V2t(x) 

+ Viu(x)Vi2T(x)-ViiT(x)Vi2i(x)] . (2.16b) 

Note that 2Q12 ^ structure similar to iQi2- This im- 
plies we could use 2Q12 instead oflQ^2- Physically, qQi2 
is related to the single particle density of states, while 
1Q12 is basically the superconducting order parameter. 



III. SADDLE POINT SOLUTIONS AND 
GAUSSIAN APPROXIMATION 

A. The saddle— point method 

We now look for a saddle-point solution of the field 
theory derived in the previous section. The saddle-point 
condition is^'^'' 
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(3.1) 



According to Eqs. (2.16), the saddle point values of both 
Q and A in singlet superconductivity-like phases have 
the structures 



r'3i2(x) = <5aia2 ^iO [^n 



^rl Qni 



r-^l2i-^) 



sp 



~t~ ^ni,n2 



(3.2a) 
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where we assume A„ = — A_„, A„ = — A_„ which is 
equivalent to a redefinition of the chemical potential, 
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and set Qn = Q-m <ln = 9-n which follows from Eqs. 
(3.2) and (2.11c). Substituting this into Eqs. (2.12) - 
(2.15), and using the saddle-point condition Eq. (3.1), 
we obtain the saddle-point equations 
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Here 



a„(k) 

^«(k) 



-(w„ " An) - 
-(*C.V.-A„)2+e^ + g2 ' 

qn 



are Green functions with = k^/2TO — \i. 
From Eqs. (3.3), it is easy to find 

k 

_2r(«)Tl^^^„(k) . 

k m 

We now define a gap function A by^^ 
9n = 9n + A = r?„A , 

with 

1 1 



(In = 



and it can be shown that 

VnU)n = i^n + i^n ■ 

We then obtain the gap equation, 

A=-2r(^)Ti5:5: 



'?nA 



(3.3a) 

(3.3b) 

(3.3c) 
(3.3d) 

(3.4a) 

(3.4b) 



(3.5a) 



(3.5b) 



(3.6a) 



(3.6b) 



(3.6c) 
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_2p(c) 
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4 + ^^ + A2 



with iV(0) = ^ the density of states per spin at the 
Fermi surface. A remarkable aspect of this gap equa- 
tion is that in this approximation the gap A and the 
critical temperature Tc are independent of the (nonmag- 
netic) disorder, and so are all thermodynamic properties 



in superconductivity. This result is known as Anderson's 
theorem. 

We next obtain the density of states. From Eq. (2.16a) 
it follows. 



7V(e^+a;) = -Re(°g„Jx) 



(3.8) 



In saddle point approximation, we have for the density 
of states 

-2 1 



— ^ 1 \ — > 

N{eF + uj) = — — > Im ^?„(k, ia;„ ^ w -h iO) 

for a; > A 



TT V 

= Nf 
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= for w < A . 



(3.9) 



For later reference we also define a matrix saddle-point 
Green function 



(3.10a) 



whose matrix elements are given by 

(Gsp)„m(k) = 5nm Gn{^) (tq «> Sq) 

-?"n(k) (ti (g) So) 



(3.10b) 



Note that the above results are the standard ones. 



B. Gaussian approximation 

Wc next set up the calculation of the Gaussian fluc- 
tuations about the saddle point discussed above. In the 
following section these results will be used to compute 
the physical correlation functions in the disordered su- 
perconducting phase. To this end, we write Q and A in 
Eqs. (2.12) - (2.15) as. 



Q = Qsp + SQ 
A = Asp SA 



(3.11a) 
(3.11b) 



and then expand to second or Gaussian order in the fluc- 
tuations 6Q and 6A. Denoting the constant saddle point 

contribution to the effective action by ^sp, and the Gaus- 
sian action by Ac, we have, to the Gaussian order, that 



(3.7) with 



A[Q,A]=Asp+Ag[Q,A] , 
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(3.12) 



Ag[Q, A] = AintlSQ] + AdisiSQ] + ^ Tr [G,p6AG,p m) 
+ j dxti (^SA(x) SQ{x)j , (3.13) 

For the quadratic part we find 
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k 1,2,3,4 r,s i,j 



V A 

rs^l2,34 



(k)i(5A)34(-k) . 

(3.14a) 



Here 



rs^l2,34 



-h] 12 



(k)=<5i3<524^?^(k)iV,","<5. 

+(5l3^2,-4'^?2(k)iV°l%;/l2 
+5i,_3 524(^}^(k)iV,l°(5,,;/ 
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+5l,_3 52,-4¥'}2(k)^;.'^i,;/l2. 



rs^l2,34 



(k);/i 



(3.14b) 



with 4x4 matrices 



iT3 
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iV" = ( 


^ -iT3 
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and 



r-'l2 



l + (5i 



+ 



(3.14c) 



(3.14d) 



where y + j = ^ro + ^ri + ^r2 — Srs, etc. and 

^n°r„(k) = ^E^"(P)^-(P + '^) ' (3.14e) 



Q and A can now be decoupled by shifting and scahng 
the A field. If we define a new field A by 

;(M)l2(k)=2j,i(A-l)i2,34(k) 

X (i(<5A)34(k) - i(<50)34(k)) ii334, (3.16) 

with being the inverse of the matrix A defined in 
Eq. (3.14b), then A and Q decouple. Integrating out the 
^A fluctuations leads to a Gaussian action completely in 
terms of SQ fluctuations, 



Ag[Q] 



-VT.J2J2J2 rimM ;i(A-l)i2,34(k) 
k 1234 rs ij 

+Ant[<5Q]+^dis[^g] , (3.17) 
It is convenient to rewrite this result as 



^g[q] = -^;r E E E E ^('^^)i2(k) ;iMi2,34(k) 

k 1234 rs ij 

xi(5Q)34(-k) , (3.18a) 



where 

;iM,234(k) = ;{(A-i)i2,34(k);Bi2i534 

-2TT^''Hij6rsSi+2,3+4: ^ 
r-^12 %^rs^l3<^24 



(3.18b) 



and (fi^^, and ip^^ deflned similarly with QQ in Eq. 
(3.14e) replaced by {—1)QJ^, {—1)J^Q, and J^J^, respec- 
tively. 

In a similar way, the term that couples SA and 6Q can 
be written 

Tr(M5Q) =4 E ^EE^('^^)i2(k) 



V 

1,2,3,4 k r,i 

x;B,2(k);((jQ)i2(-k) 



where 



r-^12(k) — r-^12 



(3.15a) 



(3.15b) 



IV. PHYSICAL CORRELATION FUNCTIONS 

A. Ultrasonic attenuation by saddle-point 
appr oximat ion 

We now use the results of the preceding sections to 
calculate transverse ultrasonic attenuation in both clean 
and disordered superconductors. As shown in Ref. 17, 
the sound attenuation coefficient has the expression 



OJ 



a{uj) = lim 
where, with D. 



■ Imx(k, iojn ^ to + iO) 



(4.1) 
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B 3 



X(k, iun) = dxdx'dydy' exp (-ik • (x - y)) ^ S{x - x')S{y - y')D^Dy 

^ (71,(72 

x;^ E(C„.,WC,-a.„,..(x')C„..(y)C.+c.„,..(y')> • 



(4.2) 



By introducing a source term of the form 
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(x), 



(4.3) 



we can obtain 



X(k = 0,iu)n) 



1 



m1(3V dh{u>n, k)9/i(ti;_„, — k) 
with the third term of the right side of the Eq. (2.12) becoming 

^3 = i Tr In (g^i -iA + D 



(4.4) 



h=0 



i Tr In (Gq ' - iiAsp + SA) + d) 



Tr In (l + DGs^ - iSAGs 



and D = S{uji — uj2 + Ct'„)/iexp (— «k • x)Dx- 

In the saddle-point approximation, we neglect the SA 
item and have 



■Trln(G7;) 



(4.5) 



B. Physical correlation functions by Gaussian 
fluctuations 



^3 = ^ Tr {DGspDGsp) + const. (4.6) 



We then obtain the ultrasonic attenuation coefficient, for 
small frequency, 



as{uj) = an 



1 + exp (/3A) 



(4.7) 



for both clean and disordered superconductors. Here a„ 
is the attenuation coefficient of the normal metal. In 
the clean metal it has 



an,clean 



with the usual conditions w < qvf < A satisfied. In the 
disordered case 



2N{0)ky'' 



^n, disordered 



(4.8b) 



where the approximation of TeA <C 1 is assumed, which is 
called the dirty limit. ^® The above result confirms the one 
of Levy's which was obtained by Boltzmann's transport 
equation. ^° It is noted that no Green function method 
has been used to obtain this result before. 

The above method can be used to obtain other phys- 
ical properties, like longitudinal electrical conductivity. 
In that case, higher-order corrections must be included 
due to the gauge invariancc problem. ^-'^ Below we show 
how to correctly obtain the conductivity by using the 
Gaussian fluctuations about the saddle point. 



1. Gaussian propagators 

We now expand our method to calculate the Gaussian 
propagators and then to obtain the number density sus- 
ceptibility, Xn, the spin density susceptibility, Xs and the 
conductivity. We find in Appendix A that the number 
density susceptibility, Xn, and the spin density suscepti- 
bility, Xs I can be expressed in terms of the Q-correlation 
functions, 



X«(k,a;„) = l|?^^ 



V 



1.2 r=0,3 



.W)i+„,i(k) 



(4.9) 



(4.8a) with x^"^ = Xn and x 



/(1,2,3) 



Xs- Here the Gaussian 



propagators in the Eq. (4.9) are given in terms of the 
inverse of the matrix M defined in Eq. (3.18b) by 



V 



"■i M 



1 

12,34 



(ki), (4.10) 



where (. . .)q denotes an average with the Gaussian ac- 
tion Ag- We see from Eqs. (4.9) and (4.10) that M"! 
determines the correlation functions within Gaussian ap- 
proximation. 

In the following section we will be interested in the 
number density susceptibility Xn- Other correlation func- 
tions can be obtained similarly by applying the technique 
introduced below. From the expression of Q in terms of 
the fermion fields, Eq. (2.8), it is easy to see that the con- 
tributions to Eq. (4.9) from r = and r = 3 are identical 
for 7^ 0. We can therefore write 

X„(k,a;„)=4T^°°M-^„^,^2+„,2(k) , (4-11) 

1>2 
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To find Ei,2 33-^i+„,i;2+n,2' ^6 rewrite M as 



= C + CDM-'^ 



(4.14) 



rs-^12,34 



with 



= ;i(c-i)i2,34(k) 



rs-^12,34 



Then we find 



(4.12) 

(4.13) 



rs*-^ 12,34 



rs^l2,34 



34 



(c-1 

It is convenient to write the inverse of the matrix M as 
an integral equation, 



J 



For further simplicity, we set T = 2TV^''\ r° = 
and \.I-^^2 = 1 fo'^ 0- Expanding Eq. (4.14) 



oo»x-i _ 00 4(0) _p 

33-'" 12,34 ~ 33^12,34 



}^Er8^-l+2,7+8 23-78,34. 



/ 



ooa#-i 



+7^ 



iM 



\ 

1,2;3,4 \ 



*^12 23'" 1,- 
+ 'f'l2 23^" -1,2;3,4 
\ ^^^12 33^*^ -l,-2;3,4 / 



(4.15) 



= irNpTe 

we have 



(4.16) 



where we have used the structures of B, Eq. (3.15b) and A^^\ Eq. (3.14b). ^ in turn obeys the integral equation 



23^^^ 12,34 



.00.(0) 
23^12,34 ^ 



,_00 f OO/i/f-l 

-'i^l2 2^7S ^1+2,7+8 23 -'^^^ 78,34 

Ml E78 ^■ 



_1_2,7+8^^M 




,„uu uu n^-^ 

-V12 23*'l,2;3,4 
V12 33^"l,-2;3,4 
^'*'l2 33^" -1,2;3,4 
2;3,4 



\ "^12 23^" 



(4.17) 



Similar results can be obtained for ^^M_-^ _2.^ 4 and 

23-"^ -l,-2;3,4- 

We can now obtain the Ei 2 i-2+ra 2 now. Ob- 

viously El,2 23^r+„,-l;2+„,2 and El,2 23^lLn,l;2+rx,2 

need to be determined first. We find in Eq. (4.17) that 

all of the 23-^ ^ form a linear equation group which can 
be solved by using Cramer's Rule. It is then easy to ob- 

tain !^§Af^;;„,i^2+.„,2 by Eq.^ (4.16). Through Eq. (4.11) 
the number density susceptibility Xn can finally be eval- 
uated explicitly. Note that this technique can then be 
generalized to obtain all elements of M~^, which in turn 
gives the Gaussian propagators or physical correlation 
functions completely. 



2. Correlation functions in the dean limit 

In this section we discuss the clean limit, or the non- 
impurity electron gas. Let us perform the clean limit, 
Te ^ 00. ^dis then vanishes. That also means A„ ^ 0, 
A„ — > and g„ = A. 

For small |k| and we obtain the number density 
susceptibility of clean superconductor. 



X„(k,a;„) = -Nf 



(4.18) 



The electrical conductivity a is determined by Xn via' 
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a-(k, w) = ie Xn(k, iujn ^ u; + iO) . 



(4.19) 



In particular, the real part of the conductivity as a func- 
tion of real frequencies has a delta-function contribution 



Re cr{uj) = - '™ 



lim e"^-^ Imxra(k, iw„ ^ u) + iO) 
fe— ►o k 



Sic,) 



mre 



m 



(4.20) 



with n = 3^ the particle number density. This coincides 
with the result already known. 

Similar procedure can be applied to obtain the spin 
density susceptibility, by noting that 



X.(k,a;„ = 0) = -rr l](^(<50)i,i(k)^(<5Q)2,2(-k) 

1,2 

1,2 



;2,2 



(4.21a) 



and 



33-'" l-|-n,l;2+ri,2 ~ / j 33^1+w.l:2-l-n.2 
1,2 1,2 

_ -Nf Un 
~ 2T n 

for ojn = 0, |k| — > 0, 



(4.21b) 
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where n = + n„ , with Us the density of superconduct- 
ing electrons, n„ the density of normal electrons, 



r 

nn = n 

J — c 



T(l + exp(^^i2j±^))2 



(4.21c) 



The result Xsik — > 0,a;„ = 0) = —Np^ is consistent 
with Yosida's.^^ 

The above result means Xs = at zero temperature. 
This is because a BCS superconductor is a perfect dia- 
magnct at T = O.^'' The non-zero part comes from the 
contribution of normal electrons at finite temperature,^* 
since some Cooper pairs are broken into normal electrons 
at T ^ 0. 



3. Correlation functions in the disordered case 

Now wc turn to the disordered case. The approxi- 
mation of TeA <C 1 is again assumed. Calculations at 
T — > show that in the limit of long wavelength and low 
frequency. 



Xn(k, Un) 



-Np- 3 



^2 + 



(4.22) 



and the real part of the conductivity as a function of real 
frequencies has also a delta-function contribution 

Re a{Lj ^ 0) = — hm e^— Imxn(k, iuin ^ co + iO) 
fe^o k 



Np A Te TT^ vj 



6{oj) . 



(4.23) 



Note that to satisfy the f-sum rule in the disordered case 
the conductivity will not vanish completely at finite fre- 
quency. Our calculation shows that at T = 0, 



Re (j{uj > 2A) = — 



dE 



-E{E - w) - 



- AV(^ - w)2 - A2 



(4.24) 



where the vertex corrections resulting from the impurity 
scattering and the interaction have been omitted. 0-„ is 
the conductivity of normal metal. This coincides with 
the result already known, too.^"'' 

Again, similar procedure can be applied to obtain the 
spin density susceptibility. We find that, at T = 



X.(fc^O,w„-0) = , 



(4.25) 



That means the spin response in the nonmagnetic disor- 
dered case is the same as that in the clean limit. This 
is consistent with Dcvcrcaux and Bclitz's argument, 
which has shown that the nonmagnetic disorder has no 
eff'ect on the spin-flip pair breaking rate. 



V. CONCLUSION 

In this paper we presented a method to study the 
transport properties of disordered s-wave superconduc- 
tor. The crucial idea is to first identify the saddle points 
of the system by using a symmetry analysis, then to study 
the fluctuations around them to obtain the physical cor- 
relation functions. The ultrasonic attenuation, number 
density susceptibility, the spin density susceptibility and 
the conductivity have been calculated in the clean super- 
conductor, as well as in the disordered superconductor. 
Other properties, like energy correlation function, can be 
similarly obtained. Furthermore, the formalism here can 
be a powerful tool to study the quantum phase transi- 
tions between normal metal and superconductor. 

Finally, we remark, that the techniques used here can 
be used to study gapless s-wave superconductors, as well 
as, for example, disordered d-wave SC relevant to the 
high Tc superconductors. 
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APPENDIX A: CORRELATION FUNCTIONS IN 
TERMS OF Q MATRICES 

The real number density susceptibility has the follow- 
ing form^ 

X^(xiti, X2t2) = -ie{ti - t2)([n(xiti), n(x2t2)]) (Al) 
where 

h = n- (n) (A2) 

with n the number density operator. It is inconvenient to 
calculate it directly. Instead, we introduce a correspond- 
ing temperature function that depends on the imaginary- 
time variables 



Xn(xiri,X2r2) = -(r.r[«(xiri)n(x2r2)]) 



(A3) 



where we have the following relation between Eqs. (Al) 
and (A3) with the Lehmann representation 



X^{]s., oj) = Xn(k, iu)n CO + iO). 



(A4) 



The time-order indication of Eq. (A3) will disappear 
in the functional integral form,^'' which is the case in the 
present paper. 

Next we notice that 
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oQnin2 — 

a 

3Qnin2 — Q ^ ] (V'ni,o-V'n2,o- ~ 4'n2,i7''Pni,a) ■ (A5b) 



By using Eqs. (A3) and (A5) we can then obtain 



1,2 r=0,3 





X "W)2+„.2(-k) 



(A6) 



Similar analysis can be applied to find the spin density 
susceptibility. With the spin density 

n,{k, ujn) = ^{ipip, w), a ip{p + k, w + w„)) (A7) 



we can obtain Eq. (4.9). 
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